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DIRECT TOMOGRAPHIC RECONSTRUCTION 



FIELD OF THE INVENTION 



The present invention relates to the field of medical imaging systems and methods and 
particularly to methods and apparatus for the acquisition of tomographic data, especially PET 
5 data and the reconstruction of three-dimensional images based on the data. 



Gamma cameras known in the art of Nuclear Medicine (NM) imaging produce 
tomographic images that are indicative of physiological activity. Such cameras receive 
radiation that is emitted by radioisotope markers or tracers, which are introduced into the body 

10 of a subject and are taken up by an organ of the body in proportion to the physiological 
activity of interest. The radiation emitted is generally received by a scintillator/detector 
system, which produces electrical signals responsive to photons of the radiation incident 
thereon. The signals are processed and back-projected, using computerized tomography 
methods known in the art, to produce a three-dimensional image indicative of localized 

15 activity within the organ. 

Positron emission tomography (PET) is a system of tomographic nuclear imaging 
which is well known. In general, this system is based on utilization of radio-isotopes which, 
during a decay, emit two photons in directly opposite directions. A ring, or rings of discrete 
detectors which surround a subject into whose body such isotopes have been introduced, 

20 detects the occurrence of such a decay by detecting two coincident gamma rays impinging at 
two detectors, where the events have an energy associated with the decay. 

Based on this coincidence detection, the position of the decay (i.e., the presence of the 
radio-isotope which decayed and caused the coincident detection) is known to be along the line 
joining the two positions at which the coincident impingements were detected. 

25 The exact origin of the event is not known and the calculation of an image of the 

distribution of the radio-isotope is, in the prior art, based on a probabilistic smearing of the 
event into sinograms associated with the line on which the event is known to have occurred. 
This smeared probability is forward projected to form (together with other detected events) 
tomographic views at each of a plurality of slice positions, generally coincident with the rings 

30 of detectors. These views are used to generate tomographic images of slices corresponding 
generally to the positions of the rings. 



BACKGROUND OF THE INVENTION 



1 
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One of the major problems with this reconstruction system is the "blurring" of events 
in the direction perpendicular to the plane of the rings. This blurring, if not reduced, results in 
an image which is unsuitable, for diagnostic purposes. 

One method of reducing the effects of axial blurring is just to reduce the "acceptance 

5 angle" for events. If events cause coincidence detection in widely spaced rings, they naturally 
cause greater axial blurring. By reducing the acceptance angle, i.e., the angle of the line 
connecting the detection points with a ring, amount of axial blurring may be reduced, at the 
price of rejecting valid detected events. 

A second method, described in U.S. Patent 5,331,553, to Muehllehner et al., axially 

10 rebins the events based on a deblurring function. This rebinning reduces the axial blurring, 
however, it increases the image noise. Substantial blurring remains and artifacts are generally 
generated. The rebinning may be performed on the views or may be performed on the three 
dimensional tomographic image. 

Another method which has been mooted is to use an expectation maximization (EM) 

15 algorithm 3-D reconstmction. Such an algorithm is described, for example, in "Maximum 
Likelihood Reconstmction for Emission Tomography" IEEE Trans Med Imag MI-1; pages 
113-122 (1982). This algorithm performs iterative expectation maximization operations on an 
equation relating the activity in the voxels with the projection data. During each iteration, all 
of the projection data is taken into account. However, this method, while theoretically 

20 possible, requires very large amounts of computation in order to reach satisfactory results, so 
that it has not been implemented commercially. 

In order to more clearly understand the operation of the invention, it is useful to review 
the prior art EM methods. 

In the EM method, a body having a variable radiation emission density is considered to 

25 be contained in a discrete cube or other discrete region and the emission density of each voxel, 
v~(x,y,z) in the region is defined as X(v). The radiation emitted by the body is detected by 
detectors surroimding the body. If two photons are simultaneously detected by the two 
detectors (which indicates that the event may be caused by a positron interaction) and a Line 
Of Flight (LOF) between the two detection coordinates intersects the region we consider this 

30 as a coincidence acquisition (more simply as an "event"). Coincidence acquisition 
reconstmction algorithms try to determine the unknown emission density distribution in the 
region given a Ust of the LOF of detected events. 
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The classic PET scanner is built of rows of detector rings whereby each detector is a 
discrete unit. All the events detected simuUaneously by the same two detectors di and dj are 
collected in a single bin Thus each bin b defines a single LOF and every event detected in 
the bin is assumed to have originated along this LOF. 

Let v=7,...,r represent voxels of the field of view and let independent variables x(v) 
with unknown mean values of emission density X(v) represent the number of unobserved 
emissions in each of the V voxels. Suppose further that if an emission occurs in voxel v it has a 
probability of p(v,b) of being detected in bin b, then p(v,b) defines a transition matrix 
(likelihood matrix) which is known. Based on the number of y'='y(b) events detected in each 
bin it is desired to estimate the total number of the imknovra distribution of events ^^'=a(v), 
v=7,..., K For each X, the observed data has the conditional probability or Hkelihood: 



where /u(b) are the mean values of the (observed) Poisson variables yf??;, that is: 



(1) 



M^)= Y.^{v)p{v,b). (2) 
\<v<V 

The maximum likelihood estimate of /I is: 

A = arg m^^P{y\X) (3) 
The optimaUty condition of the above equation , based on the log-likelihood, is: 



^ l/I IJ l^LX V I fJl V^IJ I 

,v=\,...,V (4) 



1 A yib)Mv)p(v,b) 



where P(v) is the probabihty to detect an event emitted from voxel v: 
B 

Piv)=Jlp(v,b). (5) 
b=\ 

The EM algorithm can be considered as a fixed point iterative algorithm based on (4): 

> 

'^"■'*='Za""(V)p(V,*) 

V'=l 

In 1994, a paper entitled "Accelerated Image Reconstruction Using Ordered Subsets of 
Projection Data" IEEE Trans Med Imag vol. 13, no. 4, (1994) pp. 601-609 reported the use of 



10 
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an EM algorithm using ordered sub-sets which it called OSEM. The OSEM algorithm also 
performs an EM algorithm on the equation relating the voxel activity with the projection data. 
However, in this method sub-sets of the projection data, rather than the full set of projection 
data, are taken into account in each iteration with a different sub-set being taken into account 
for each iteration. If the projection data is divided into N sub-sets (which together form the 
complete projection data set) and each iteration is performed using only one of the sub-sets, 
then, if N iterations of this type are performed, the overall results will be comparable to those 
achieved when N iterations are performed taking into account all of the projection data. 

This surprising result allows the practical use of EM reconstruction algorithms (in the 
OSEM form) in PET systems using rings of detectors. A description of the OSEM method as 
in the above referenced paper by Hudson, et al, follows: 

hctyf be the number of photon emissions recorded in the i^^ projection bin and let Yff 
be the set of parallel projections {y^, y/+ j,,..} that view the object at angle 0, orthogonal to the 
tomographic axis. The projection data are grouped into n subsets Sj, iS'2,...,5^. If there are a 
15 total of P projection bins, the elements of each subset, in no particular order, are: 

" ^(n-l)(P/rt) + l 'y(n'J)(P/n)'h2''"yp^' 

The subsets normally consist of projection views separated by some fixed angle about 
the object. For example, each subset might consist of two sets of parallel projection views 
spaced 90^ apart, e.g., Sj^{Yo.Y^2l S2={Y^4.Ys^4J, and so on. 
20 In the OSEM algorithm, the log-likelihood objective function for each of the subsets is 

increased with each iteration, using the results of the previous iteration as the starting point. 
Therefore, the EM iterations become: 

^ TKa ^ ~V ,v=l,...,F . (8) 

V*=l 

where is the estimated nimiber of emissions from v after the introducing the I^^ subset of 
25 projections. 



4 
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The EM procedure is repeated until all n subsets have been exhausted. Such a cycle is 
considered a single iteration of OSEM. The cycle can be repeated iteratively until a 
satisfactory reconstruction is obtained. 

However, ring type PET systems are not of general utility. Their use of separate 
5 detectors for each detection pixel results in a system which can have high sensitivity, but not 
high resolution. For this reason and, to a lesser extent, because of its geometry, such a system 
is useful only for PET and cannot be used for other NM applications such as for acquiring 
planar images and for SPECT. 

Neither the EM nor OSEM algorithms are easily applied to systems of planar detectors. 

10 Such application would be very desirable since this would allow for the use of standard 
rotating dual head gamma cameras for PET as well as for SPECT. However, this is not 
practical. If an attempt is made to use a pair of opposed planar detectors (of the type normally 
used for planar or SPECT imaging), rotating about the subject, to acquire data for constructing 
PET images using EM (or OSEM) techniques, the amount of data which is acquired is 

15 reduced, as compared to the ring type PET system, making the planar system impractical. For 
example, if two detectors having dimensions of 540x400 mm and a rotation radius of 350 mm 
are used, with a bin resolution of 2.5x2.5 mm, there are about 3x1 0^ bins. If the system is 
capable of an acquisition rate (for coincident photon events) of 1000 events per second, in a 
typical study of 30 minutes approximately 1.8 million events can be collected. The number of 

20 bins is more than two orders of magnitude larger than the nimiber of events, which means that 
most bins will remain empty and almost no bins will have multiple events. Thus, for such a 
device, the normal binning procedure according to a fixed set of a very large number of bins is 
impractical and results in excessive calculation. Also, with such sparse data the results can be 
expected to be noisy. 

25 It is understood that for the classical ring type PET imager, the number of bins is much 

smaller and the absorption of the photons is more efficient (since the crystal is thicker) such 
that, OSEM reconstruction becomes practical. The prior art does not teach any way to 
combine the advantages of the OSEM system with the higher intrinsic resolution of planar 
detectors. 

30 



5 
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SUMMARY OF THE INVENTION 

Accordingly, it is an object of some aspects of the present invention to combine the 
advantages of the EM system, and in particular the advantages of the OSEM system with 
advantages of detection utilizing planar detectors. 
5 It is an object of some aspects of the present invention to provide a method of 

Tomographic images which generates a three dimensional image without generating 
intermediate projection data or sinograms. While the method is especially applicable to PET 
using planar detectors, it is applicable to all types of three dimensional nuclear tomography, 
such as SPECT and to other types of tomography such as X-Ray CT. It is also applicable to 
10 PET using ring detectors, especially those having a large number of rings. 

It is an object of some aspects of the present invention to provide a method of three 
dimensional reconstruction of images by which a three dimensional image is generated 
without first providing a stack of two-dimensional data sets representative of individual slice 
images. 

15 It is an object of some aspects of the present invention to provide a method of three 

dimensional reconstruction of PET images by which a three dimensional image is generated in 
which an Expectation Maximization (EM) methodology is used for data which is acquired 
using area detectors. 

It is an object of some aspects of the present invention to provide a method of three 
20 dimensional reconstruction of PET images using a methodology similar to that of the Ordered 
Set Expectation Maximization (OSEM) method. 

However, since the sets in the present invention do not need to be ordered and may 
even consist of single events, we prefer to refer to the methodology of the present invention as 
an Event Driven Expectation Maximization (EDEM) methodology. In particular, in some 
25 aspects of the invention, the sub-sets are not ordered according to their geometric properties. 

It is an object of some aspects of the present invention, to provide a method of three 
dimensional reconstruction of images using an EDEM methodology in which the elements of 
the sub-sets each have only a single event. Furthermore, in some aspects of the invention, the 
elements of the sub-sets include only those elements for which an event has occurred. 
30 It is an object of some aspects of the present invention to provide a method of 

tomographic reconstruction in which three dimensional images can be reconstructed using data 
having less than 180 degrees of view. As opposed to normal back projection methods, such a 



• 

wo 98/47103 



PCT/IL97/00128 



reconstruction has few if any artifacts. It is also possible to begin the reconstruction of the 
images prior to the acquisition of a "full set" of data, starting even from the first event. This 
allows for a true evolving image, i.e., a reconstructed image which is updated as data is 
acquired. This type of reconstruction is applicable to many types of gamma ray and x-ray 
5 tomography. 

Furthermore, since the method is tolerant of incomplete data (in the slice back- 
projection tomography sense), it is possible to apply the EDEM methodology to cone beam 
tomographic imaging. 

As utilized herein in the specification and claims, a spatially continuous detector is 

10 defined as a detector in which interactions are detected on a spatially continuous scintillator or 
the like and in which the position of the interaction with the scintillator are determined by 
utilizing outputs of a plurality sensors. Such detectors give generally continuous values of 
position. One example of a spatially continuous detector is an Anger Camera. 

As utilized herein the term substantially planar detector is a detector the detection 

15 surface of which is substantially flat. Examples of such detectors are conventional Anger 
cameras and solid state mosaic cameras, as known in the art. 

In a preferred embodiment of the invention at least two area type detectors are placed 
on opposite sides of a subject in whom a distribution of a radiopharmaceutical is to be imaged. 
These two detectors are rotated, generally in a continuous manner or optionally in a stepwise 

20 manner, to acquire images from a plurality of directions about the subject. As is usual when 
area type detectors are used, an Anger type methodology is utilized for determining the 
position of an interaction on the detector. When such a detector is used, the efficiency of 
detection of photons is much lower than for the thick detectors used in the rings of individual 
detectors traditionally used in PET. Also since, in general, only two detectors are used, many 

25 of the events which could in principle be used to produce the image do not reach the detectors. 
In view of these limitations, the data which is detected is too sparse to be used with ordinary 
EM or even OSEM methods. 

In a preferred embodiment of the invention, the individual positron events (i.e., pairs of 
concurrently detected interactions on opposite detectors) each forms a single element of the 

30 data which is to be subjected to the EM algorithm. In contrast to the prior art, in which events 
are binned in bins which generally include a large number of events and in which events 
having somewhat different angles and positions of impact on the detectors are grouped, the 
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present invention, in its most preferred form, has only one event in each element of the 
reconstruction, such that each reconstruction element is more sharply defined in space than in 
the prior art. 

In particular, the prior art bins events geometrically and then (in the OSEM method) 
utilizes sets of such geometrical groupings of the events to reconstruct the three dimensional 
image. In preferred embodiments of the present invention, the events are not binned 
geometrically and the sub-sets which are used in an OSEM type method need not be 
geometrically related or even related in any way. 

In order to avoid confusion with the OSEM method, the present method, which 
utilizes sub-sets which are not necessarily ordered at all, the general methodology of the 
present invention is referred to herein by the term Event Driven Expectation Maximization 
(EDEM) method. 

In a preferred embodiment of the invention, the elements of a sub-set bear a time 
relationship to each other and sub-sets of data are applied to the EDEM algorithm according to 
when they were acquired. Thus, in contrast to many earlier methods of reconstruction, in 
which all the data is acquired before the image is reconstructed, in some aspects of the present 
invention, the image reconstruction is started, using data which is acquired during an initial 
period. This initially acquired data is used to construct a first estimate of the image while 
further data is acquired. This fiirther data is acquired during a second period and this further 
data is applied to the first estimate to determine a second estimate. Data acquired during the 
period in which second estimate is determined is then used to generate a third estimate and so 
forth. 

The estimates of the image may be displayed and/or analyzed so that the operator (or 
an automatic method) may follow the progress of the acquisition and end the acquisition when 
the image quality has reached a desired level or does not improve beyond a certain point. 

In a preferred embodiment of the invention, as each event is detected (i.e., when two 
simultaneous interactions are detected on the opposing detectors) the position of the detected 
interactions is determined to a precision which is consistent with the precision available from 
the Anger camera which forms the basis for the detection of the event. This precision is much 
higher than that available from the ring type detectors normally used for PET which provide a 
single crystal and associated light detector for each pixel. 



8 
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A line connecting the two determined points at which the interactions have occurred is 
determined and the event is weighted (adjusted) for the probability that an event occurring 
along the line would be detected, for attenuation and for the effect of the angle at which the 
event is detected on the density of events which are detected. The adjusted weight of the event 
5 is then distributed along the line into voxels in the region to be reconstructed, where, 
preferably, the amounts which are distributed are proportional to the length of the line which is 
contained within the individual voxels. Since the event is well locahzed on the detectors and 
since each event is separately considered, the events can be more precisely and simply 
distributed among the voxels than in the prior art where a bin contains events with different 

10 LOFs bunched around a central LOF, 

After a certain time, which may be a preset time or which may be determined by 
acquisition of a certain number of events, but is most preferably determined by a preset 
rotation of the detectors about the subject, the events are applied to an initial distribution of 
events using an expectation maximization algorithm as is well known in the art. While the 

15 choice of an initial distribution is arbitrary, a uniform distribution is generally used, although 
other distributions, determined by the location of the organs being imaged, may be preferred 
under some circimistances. 

The resulting estimated distribution is the first iteration of an EDEM methodology 
which is applied to the data. Subsequent sub-sets of individual events are then applied, 

20 preferably sequentially, to the previous iteration until the image quality as judged by the 
operator, or by a quality determination algorithm in the camera, reaches a desired level. 
Alternatively, the acquisition may be stopped using the usual preset time, count number and 
minimum count rate criteria. 

In a preferred embodiment of the invention a smoothing, median or other noise 

25 reducing filter is applied to the data generated by a given iteration prior to the application of a 
subsequent iteration. When a filter of this type is so applied the effect on the resolution of the 
final image is small, since subsequent iterations "correct" the smoothing, while the noise level 
of the image is reduced. Alternatively, a band enhancing filter, such as a METZ filter which 
also reduces noise can be used. 

30 As indicated above, if the LOF can be determined, as for example by consideration of 

the position of the event on the detector and the collimator angle for SPECT and by 
consideration of the focal spot of the x-ray tube and the point of detection of the x-ray photons 
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(for X-ray CT tomography), the OSEM algorithm of the present invention is appHcable. 
Furthermore, since the method of the present invention does not require the construction of 
sinograms and/or slice data, it allows for direct reconstruction of three dimensional images 
from cone beam based data. Since the method is tolerant of incomplete data sets, cone beam 
5 data (generated from either single or multifocal sources) can be used, and used in an efficient 
manner to generate three dimensional data. Furthermore, fan beam data (either single or 
multifocal sources) can be used without binning and without the need for forming sinograms. 

There is thus provided, in accordance with a preferred embodiment of the invention, a 
method of reconstructing tomography images comprising: 
10 acquiring data on individual radiation events; 

distributing a weight of the individual radiation events along a line of flight associated 
with the event determined from the acquired data; and 

iteratively reconstructing the image based on the individual reprojected data. 
Preferably, the weights are distributed in voxels along the line of flight and wherein the 
15 weight of a particular event is distributed based on the probability that an event occurred in 
particular voxels. 

ha a preferred embodiment of the invention, the line of flight of an event is determined 
based on the position at which the event was detected on a detector and the acceptance 
direction of a collimator through which the detector receives radiation associated with the 
20 events. 

In an alternative preferred embodiment of the invention, the line of flight of an event is 
determined by the position on a detector on which the event is detected and the location of the 
source of radiation associated with the event. 

In an alternative embodiment of the invention the line of flight associated with an event 
25 is detennined by detection of two coincident photons. 

Preferably, reconstructing the image comprises applying an iterative expectation 
maximization (EM) method on the data in sub-sets. In a preferred embodiment of the 
invention, the individual events form the separate sub-sets. Preferably, the sub-sets are formed 
based on the time of acquisition of events. Alternatively, the sub-sets are formed from 
30 unrelated events. 

There is further provided, in accordance with a preferred embodiment of the invention, 
a method of reconstructing tomography images comprising: 

acquiring data on individual radiation events; and 
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applying an iterative expectation maximization (EM) method on the data in sub-sets 
which are formed based on the time of acquisition of the events. 

Preferably, the subsets consist of data having less than a 180 degree view angle. 
Preferably, iterations of the EM method are perfomied prior to the acquisition of data having a 
5 180 degree angle of view. Preferably, iterations are commenced on receipt of the first detected 
event. 

In a preferred embodiment of the invention, the includes displaying an evolving image 
based on successive iterations iterative method on a display device. Preferably the method 
includes determining if a study should be terminated based on the image quality of an image 
1 0 after an iteration. 

In a preferred embodiment of the invention, intermediate images are filtered with a 
smoothing fiher between iterations of the EM method. The images can be filtered using a 
smoothing or other noise reducing filter. 

In a preferred embodiment of the invention, data is reused in subsequent iterations of 
15 the EM algorithm. 

There is further provided, in accordance with a preferred embodiment of the invention, 
a method of reconstructing tomography images comprising: 

acquiring data on individual radiation events; and 

iteratively reconstructing a three-dimensional image fi-om said individual radiation 
20 events without producing a stack of two dimensional data sets. 

There is fiurther provided, in accordance with a preferred embodiment of the invention, 
a method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

iteratively reconstructing a three-dimensional image without producing individual 
25 sinograms for sUces of the three dimensional image. 

There is further provided, in accordance with a preferred embodiment of the invention, 
a method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

iteratively reconstmcting a three-dimensional image utilizing the individual radiation 
30 events without spatially binning of the events. 

Preferably, reconstructing the image comprises utilizing an expectation maximization 
(EM) method acting on individual unbinned events. 

In one preferred embodiment of the invention, the radiation events are nuclear emission 
events and the images are emission tomography images. 
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In an alternative preferred embodiment of the invention, the radiation events are 
positron annihilation events and the images are PET images. 

In another alternative embodiment of the invention the radiation events are represented 
by photons which have passed through a subject and the images are transmission tomography 
5 images. In one variation of this preferred embodiment the radiation events are nuclear 
disintegrations and the images are nuclear transmission tomographic images. In another 
preferred embodiment of the invention the radiation events are X-rays and the images are X- 
ray CT images. 

In a preferred embodiment of the invention, the lines of flight associated with the 
10 radiation events form a fan beam. 

In an alternative preferred embodiment of the invention the lines of flight associated 
with the events form a cone beam. 

There is further provided, in accordance with a preferred embodiment of the invention, 
a method of reconstructing positron emission tomography (PET) images comprising: 
15 acquiring data on individual positron emission tomography events utilizing a plurality 

of spatially continuous area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method. 

Preferably, the spatially continuous detectors are substantially planar detectors. 

There is further provided, in accordance with a preferred embodiment of the invention, 
20 a method of reconstructing positron emission tomography (PET) images comprising: 

acquiring data on individual positron emission tomography events utilizing a plurality 
of substantially planar area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method. 

Preferably, the consists of two such detectors. 
25 The invention will be more clearly xmderstood from the following description of the 

preferred embodiments thereof, taken in conjunction with the following figxu-es in which: 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a simplified perspective representation of a portion of a PET imaging and 

reconstruction system, in accordance with a preferred embodiment of the invention; 

30 Fig. 2 is a simplified cross-sectional representation of a portion of the system of Fig. 1 ; 

Fig. 3 is a simplified block diagram of a front end useful in the system of Figs, 1 and 2; 

and 

Fig. 4 is a simplified representation of the position of an event used for PET 
reconstruction, in accordance with a preferred embodiment of the invention. 

12 
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PREFERRED EMBODIMENTS OF THE INVENTION 



Figs. 1 and 2 are simplified perspective and cross-sectional representations of a PET 
imaging and reconstruction system 10 in accordance with a preferred embodiment of the 



invention. 



10 



15 



20 



25 



PET system 10, in the preferred embodiment shovra in Figs. 1 and 2 comprises a pair 
of area detectors 12 and 14 placed on opposite sides of a subject 13 who has been previously 
injected with a material which decays producing a pair of photons which are ejected in 
opposite directions. Area detectors 12, 14 can be any type of area detectors known in the art. 
However, in a preferred embodiment of the invention, the area detectors comprise a large 
scintillator crystal which emits light when a gamma ray photon is absorbed by the crystal. A 
series of photomultipliers, preferably arranged in a hexagonal configuration are attached to and 
view the crystal and produce electrical signals which are proportional to the amount of light 
which reaches the photomultiplier. The position on the crystal at which the ray or photon is 
absorbed is computed from the signals by a process known as Anger arithmetic or by other 
process as known in the art which compute a continuum of positions of events on the crystal. 
Alternatively, but generally less optimally, the area detector may comprise a mosaic of small 
scintillator crystals each of which has a light (or other interaction) detector associated with it. 
One preferred area detector of this type is described in PCT Application PCT/IL96/00164, 
filed November 24, 1996. 

When an interaction with a scintillator crystal is detected by one of the area detectors, 
the signal which is generated is passed via one of lines 16, 18 and a front end 20 to a computer 
22. While computer 22 is shown as a PC, generally, a more powerful computer is required to 
perform the calculations described below in real time or to complete them in a reasonable time 
after data acquisition is completed. 

The pair of area detectors 12, 14 rotates about the subject as shown by arrow 23 such 
that sufficient information is acquired from all directions to construct a tomographic image. 

Fig. 3 shows a block diagram of front end electronics 20. Front end electronics 20 
receives a signal from one or more of lines 16 and 18. A coincidence detector 24 determines if 
signals having a specified energy that is characteristic of the positron emission, are received 
coincidentally on lines 16 and 18. If the signals meet the requisite energy and coincidence 
requirements, the signals are passed to position calculation circuitry 26 of any of the types 
which are well known in the art, to determine the position of the interaction. In a preferred 
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embodiment of the invention, the position is determined to a high accuracy such that each 
detected positron emission is characterized by a pair of locations on the scintillator crystals 
which is almost always different from that of any other detected emission. As will be 
described below, the chances that two emission events have the same pair of coordinates (and 
5 LOF) is small. If two should happen to have the same coordinates they can be treated 
separately, as though their coordinates were different. Alternatively, the positions of the 
interactions may be found before coincidence is determined. Any of the operations described 
herein can be performed in software (after the signals are digitized) or in hardware as is known 
in the art. 

10 It is a characteristic of some preferred embodiments of the invention that the events are 

not binned to have a common pair of locations and that each event is considered separately 
from other events. 

In a preferred embodiment of the invention, the position of interactions are considered 
with respect to a confining cylinder 28, the cross-section of which is shown in Fig. 2, which 

15 has a length generally the same as the extent of detectors 12 and 14 in the direction of the 
height of subject 13. In general, the calculations required for the method described below (and 
for other methods of PET image reconstruction, including the standard EM and OSEM 
methods, described above) are substantially simplified if each event is considered to have been 
detected on confining cylinder 28 rather than on detectors 12 and 14. This is especially true of 

20 the calculation of P(v). Utilizing the confining cylinder considerably simplifies the 
calculations, however, the method of the invention can be applied utilizing methods of the 
prior art for calculating P(v), 

Fig. 4 shows the geometric conventions utilized in the following description of the 
preferred embodiment of the invention. Confining cylinder 28 is denoted in the following 

25 discussion as S, where: 

S:x^ ^y^ ^ R^-H<z<+H (9) 
where the cylinder has a height 2H (generally equal to the length of the detector along the 
length of the subject) and a radius R, where R is radius of the circumscribing cylinder 
described above. 

30 In the following calculation, the events are considered to be collected on a confining 

cylinder of perfect collection efficiency and not on the detectors. This transformation is 
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performed by multiplying every detected event by a profile correction factor, described below, 
and/or by other correction factors. 

Planar detectors have a less than perfect physical and spatial gamma ray capture 
efficiency. Consequently, a corrective event weight is required. Consider a set of planar 
detectors, covered by multi layered graded absorbers to shield it from events degraded by 
Compton scatter inside the imaged object (low energy filter). The capture efficiency of a 
coincidence gamma pair along an LOF located at a distance x fi*om the center of rotation and 
inclined at an azimuthal angles and rotational angle y/ with respect to the detector surface is: 

-i2 



«=\nax(^) 



where is the path length traversed by the LOF in layer i of the graded absorber, and ju^ is its 
gairuna attenuation coefficient at the energy of the photons, generally 511 KeV for PET. 
Similarly, T((p, y/) is the pass length traversed by the LOF in the detector and /j is the detector's 
gamma attenuation coefficient at 511 KeV. ^oinax(^) determined by the detectors 
dimensions. Consequently, in order to compensate for the partial capture efficiency of the 
15 detection system, each event is assigned a weight of , This is called profile correction. 

Since events are considered on an event by event basis, each event can be corrected, for 
example, for attenuation, scatter and other factors can be included on an event by event basis. 
For example, scatter can be included by reprojecting events along a cylinder having the line of 
flight of the event as its axis (optionally, with varying weight of distribution in a transverse 
20 direction depending on the distance from the cylinder axis) rather than along the line of flight 
itself. Similarly attenuation information or assumptions can be used to adjust both P(v) and 

Consider, without loss of generality, a point p=(r,0,h) inside the cylinder having a 
distance r firom the cylinder axis and a shift h in axial direction fi-om the center of the cylinder. 
25 Let the line of flight of a positron emission originating at point p be defined by (r, 0, K (p. y/). 
where: 
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0<r<R 
0<h<H 

-K<(p<n: (10) 

Then using a straightforward geometric calculation, the probabiHty that an event at 
coordinates (rM), for any angle, will be detected are: 
n 



p{r,h) = - J-==£^= ,0<r<i?,0<A<// (11) 

^ovi + ^ (^) 



where A(<p) — 



r -cos^H- 

-r-coscp + ^fn h r 

— , cp^<^<n, -<- (12) 

-r - Qostp^^lD h r 

where D-R^-r^sin^tp. 

The integral is conveniently computed using the Simpson integration formula. The 
interval (0,^) is divided, for example, into 100 subintervals. This typically provides an 
accuracy of 10'^. Note that in the third case of equation 12, cpo should be one of the discrete 
integration points. (Pq is given by: 

(p = arccos . 0< < - 

This probability of detection, which is computed using equations 11 and 12 can be 
used for the term P( v) in the EM and OSEM methods of determining the distribution of events 
described above as well as in the preferred embodiments of the present equation. 

In a preferred embodiment of the present invention, a EDEM method, which is similar 
to the OSEM method described above, for ordered sets, is used. Unlike the sub-sets of the 
standard OSEM method, the sub-sets of the present invention are preferably built according to 
a listing of unassociated events or events whose only association is that they occurred during a 
given time period. In particular, the present invention does not, in its most preferred 
embodiments, use projections or sinograms. 
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Let L represent the listing of events. L is split into m sub-sets of (preferably sequential) 



events 5/, S2,'^-S^ such that the sub-sets are roughly equal in size. For example, an 
appropriate partition would be to assign to each sub-set the events detected in a half or full 
rotation of the detectors about the subject. This sub-set of events is, in principle, sufficient to 
generate a tomographic image. In this partition, m would be the number of rotations (or half 
rotations) completed during the acquisition process. Assuming that each rotation takes 3 
minutes, there would be about 15 or 30 sub-sets in a typical study. 

However, it should be understood that the present method can be applied utilizing sub- 
sets which do not include what would normally be considered a "complete set" of data, namely 
data having views over at least 180 degrees. In fact, the data at a single position of the 
detectors can form a sub-set. Furthermore, the method can be appUed where the subsets are 
formed utilizing a single event. Thus, for the present method, the evolution of the image can 
begin starting with the capture of the very first event. 

The EM formula used to apply the sub-sets in successive iterations is, for example: 



where N is the total number of events in S, Nj, is the number of events in the j^^ sub-set, k is 
the number of the previous outer iteration and w(b) is a profile correction factor assigned to 
event b. Note that since each event is considered separately, a different profile correction is 
calculated and assigned for each event. It should be clear that the last inner iteration gives the 
starting point for the next outer iteration. It should be imderstood that while this form is 

convenient for calculation the term Aj may be removed fi-om the summation, since it is a 
constant for the summation. 

A major advantage of the present invention is that iterations are performed using 
elements of the possible data sets which actually have events associated with them. Thus, in 
this sense, the method is computationally optimal, since no null data points are considered. 

In a preferred embodiment of the invention a smoothing, median or other noise 
reducing filter is applied to the data generated by a given iteration prior to the application of a 
subsequent iteration. When a filter of this type is so applied the effect on the resolution of the 




xUv)p{v,b)iv(b) 



(13) 
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final image is small, since subsequent iterations "correct" the smoothing, while the noise level 
of the image is reduced. 

A major problem in utilizing any of the three dimensional EM method of the prior art 
in the context of a dual head rotating scanner is the need to bin the data into sinograms or 

5 projections. To apply the method to the acquired data would require either using a low 
resolution (compared to that available jfrom the camera) which would result in lower quality 
images or to use a huge number of sinogram bins (most of which are empty, in any event) 
which will make the method very inefficient. The present invention provides a method which 
combines the advantages of the OSEM approach without the need for excessive computation. 

10 It should be understood that the methods of the present invention are also applicable to 
standard ring type PET systems and may be very useful when a large axial aperture and a large 
number of bins are available. In this case, the events are grouped into time-based subsets or 
individual events form the sub-sets as in the case where a large area detector is used. 



The preferred method of the present invention is thus seen to include two parts. First, 



15 P(v)=p(r,h) must be calculated. Then the iterations described above must be performed. The 
main computational work to be performed is the ray tracing of the events, which determines 
the distribution of the probability of the event into the voxels which the ray intersects. This 
distribution is denoted by p(v,b) or p(v\b) in the above equations, it being understood that the 
distribution takes into accoimt weighing factors for each event As indicated in equation 13, 

20 these probabilities are multiplied with the current estimate to create the new estimate. A layout 
of an algorithm useful for implementation of the method is: 



25 Iteration: 

for every outer iteration k=J,2,3., do 
for subsets S(j)J=L,,m do 
N0)=0 

for every event b in S(]) do 
30 s(b)^0, N0)'=N0)-^1 



Initialize: 



1 . create the probability vector P(v), 

2. fix the initial estimate of ^ and N, 



for every voxel v intersected by b do 
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compute intersection length p(v,b) (ray tracing) 
compute t(v) = A,^ (v)p{v,b)Xb) = s(b)-ht(v) 



enddo 



for every voxel v intersected by b do 



5 




enddo 



enddo 



calculate the next estimate: /l*. . (v) = 



NU) P{v) 



x(v) 



enddo 



10 



set A/\ — ^TH 



enddo 



As indicated above, the most computational intensive portion of the iterative procedure 
is the computation of t(v), the weight of the event to be added to voxel v. An efficient method 
for finding this factor, follows: 
15 Rav Tracing 

Ray tracing starts fi-om the given pair of points which characterize the event, (xl*, yl*, 
zl*) and (x2*, y2*, z2*) which represent the intersection of the LOF of the event onto the 
surface of the field of view (FOV). The coordinates are measured in voxels and should be 
positive inside the field of view. The following is a parametric expression for the ray: 
20 x=xi*+aT, y=yi*4-bT, zi*-+-cT 



(19) 



where 



0<T<L, a-(x2*-xi*)/L, b=(y2*-yi*)/L, c=(zi*-z2*)/L (20) 



where L = ^(^2 ■^(y2 *-3^i*)^ +(^2 • 



25 



The aim of the ray tracing is to compute the length of intersection between the ray and 
the voxel. Let int(w) be the integer part of a real number u and let: 



7=int(y)+l b>0; y=int(y) b<0; z=int(z)+l c>0; z=int(z) c<0; ^=int(x)+l 
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This last term is always true since a>0 can always be obtained by approximately 



determining which of the intersection points will have an index 1 and which will have an index 
2. 

The ray tracing algorithm treats the voxel space as a grid of integer numbers. To find 
when a ray hits a grid point it is sufficient to calculate the first time x, y or z becomes an 
integer. Using the above notations the ray tracing algorithm can be described as follows: 

set T=0, let sbe a small number. 
While T<L, do 

x—xj'^+aT; y—yj'^-^bT; z—zj^-^cT. 

i=int(xj; j=int(y); k==int(z) 



The case where a=0, b=0, c=0 is avoided by adding a small nimiber to the values. 

The calculation formula for (;c,j;) depends only on the values of (b,c). These are 
determined before the main loop has started. By separating the program into four different 
loops, all conditional statements are eliminated. This makes the algorithm highly efficient. 

It should be understood that in view of the fact that intermediate images are readily 
available the reconstruction region over which the weight probability of events is distributed 
can be adjusted as the image develops, resulting in more efficient and accurate image 
reconstruction. Thus, as the image evolves, the weights of the events may be distributed only 
over regions which actually contain radiation sources. 

While the invention has been described with respect to a preferred embodiment thereof 
utilizing two area detectors having a single crystal scintillator, various aspects of the invention 
are also applicable to other types of gamma cameras, for example, cameras with multiple 
crystal area type detectors and with ring type gamma cameras of the type normally used for 
PET. While, with these gamma cameras, the resolution will be lower than for single crystal 




T=T+t 



end while 
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Anger type cameras, the advantages of on-line reconstruction and viewing of images based on 
partial data are available, using the method of the invention, for these types of cameras. 

Other variations on the preferred embodiments of the invention will occur to persons of 
the art. The present invention is thus not to be construed as being limited by the preferred 
embodiment which is presented as a non-limiting example only, but rather the invention is 
defined by the claims in which: 
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CLAIMS 

1 . A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; 

distributing a weight of the individual radiation events along a line of flight associated 
5 with the event determined from the acquired data; and 

iteratively reconstructing the image based on the individually reprojected data. 

2. A method according to claim 1 wherein the weights are distributed in voxels along the 
line of flight and wherein the weight of a particular event is distributed based on the 

10 probability that an event occurred in particular voxels. 

3. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined based on the position at which the event was detected on a detector and the 
acceptance direction of a collimator through which the detector receives radiation associated 

15 with the events. 

4. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined by the position on a detector on which the event is detected and the location of the 
source of radiation associated with the event. 

20 

5. A method according to claim 1 or claim 2 wherein the line of flight associated with an 
event is determined by detection of two coincident photons. 

6. A method according to any of the preceding claims wherein iteratively reconstructing 
25 the image comprises applying an iterative expectation maximization (EM) method on the data 

in sub-sets. 

7. A method according to claim 6 wherein the individual events form the separate sub- 
sets. 

30 

8. A method according to claim 6 or claim 7 wherein the sub-sets are formed based on the 
time of acquisition of events. 

9. A method according to claim 6 wherein the sub-sets are formed from unrelated events. 

35 
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10. A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

applying an iterative expectation maximization (EM) method on the data in sub-sets 
which are formed based on the time of acquisition of the events. 



11. A method according to any of claims 6-10 wherein the subsets consist of data having 
less than a 180 degree view angle. 

12. A method according to any of claims 6-11 wherein iterations of the EM method are 
10 perforaied prior to the acquisition of data having a 1 80 degree angle of view. 

13. A method according to any of claims 6-12 wherein iterations are commenced on 
receipt of the first detected event. 

15 14. A method according to any of claims 6-13 comprising displaying an evolving image 
based on successive iterations iterative method on a display device. 

15. A method according to any of claims 6-14 and including determining if a study should 
be terminated based on the image quality of an image after an iteration. 

20 

16. A method according to any of claims 6-15 wherein intermediate images are filtered 
with a smoothing filter between iterations of the EM method, 

17. A method according to any of claims 6-15 wherein intermediate images are filtered 
25 with a noise reducing filter between iterations of the EM method. 

18. A method according to any of claims 6-17 wherein data is reused in subsequent 
iterations of the EM algorithm. 

30 19. A method according to any of the preceding claims wherein the image is a three 
dimensional image. 

20. A method of reconstructing tomography images comprising: 



5 



acquiring data on individual radiation events; and 
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iteratively reconstructing a three-dimensional image from said individual radiation 
events without producing a stack of two dimensional data sets. 

21. A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

iteratively reconstructing a three-dimensional image without producing individual 
sinograms for slices of the three dimensional image. 

22. A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

iteratively reconstructing a three-dimensional image utilizing the individual radiation 
events without spatially binning of the events. 

23. A method according to any of claims 20-22 wherein reconstructing the image 
comprises utilizing an expectation maximization (EM) method acting on individual unbinned 
events. 

24. A method according to any of the preceding claims wherein the radiation events are 
nuclear emission events and the images are emission tomography images, 

25. A method according to any claims 1-24 wherein the radiation events are positron decay 
events and wherein the images are PET images. 

26. A method according to any of claims 1-24 wherein the radiation events are represented 
by photons which have passed through a subject and wherein the images are transmission 
tomography images. 

27. A method according to claim 26 wherein the radiation events are nuclear 
disintegrations and wherein the images are nuclear transmission tomographic images. 

28. A method according to claim 26 wherein the radiation events are X-rays and wherein 
the images are X-ray CT images. 

29. A method according to any of the preceding claims wherein the line of flight associated 
with the radiation events form a fan beam. 
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30. A method according to any of claims 1-28 wherein the lines of flight associated with 
the events form a cone beam. 

31 . A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on individual positron emission tomography events utilizing a plurality 

of spatially continuous area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method. 

32. A method according to claim 30 wherein the spatially continuous detectors are 
substantially planar detectors. 

33. A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on individual positron emission tomography events utilizing a plurality 

of substantially planar area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method. 

34. A method according to any of claims 31-33 wherein the plurality of detectors consists 
of two such detectors. 

35. A method according to any of claims 31-34 wherein the images are three dimensional 
images. 
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[received by the International Bureau on 12 Oune 1998 (12.06.98); 
original claims 1, 31 and 33 amended ; 
remaining claims unchanged (2 pages)] 

1 . A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; 

separately distributing a weight of each of the individual radiation events along a line 
5 of flight associated with the event determined from the acquired data; and 
iteratively reconstructing the image based on the distributed events. 

2. A method according to claim 1 wherein the weights are distributed in voxels along the 
line of flight and wherein the weight of a particular event is distributed based on the 

10 probability that an event occurred in particular voxels. 

3. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined based on the position at which the event was detected on a detector and the 
acceptance direction of a collimator through which the detector receives radiation associated 

15 with the events. 

4. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined by the position on a detector on which the event is detected and the location of the 
source of radiation associated with the event. 

20 

5. A method according to claim 1 or claim 2 wherein the line of flight associated with an 
event is determined by detection of two coincident photons. 

6. A method according to any of the preceding claims wherein iteratively reconstructing 
25 the image comprises applying an iterative expectation maximization (EM) method on the data 

in sub-sets. 

7. A method according to claim 6 wherein the individual events form the separate sub- 
sets. 

30 

8. A method according to claim 6 or claim 7 wherein the sub-sets are formed based on the 
time of acquisition of events. 

9. A method according to claim 6 wherein the sub-sets are formed from unrelated events. 

35 
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30. A method according to any of claims 1-28 wherein the lines of flight associated with 
the events form a cone beam. 

5 31. A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on individual positron emission tomography events utilizing a plurality 

of spatially continuous area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method acting on 

individual xmbiimed events. 

10 

32. A method according to claim 30 wherein the spatially continuous detectors are 
substantially planar detectors. 

33. A method of reconstmcting positron emission tomography (PET) images comprising: 
15 acquiring data on individual positron emission tomography events utilizing a plurality 

of substantially planar area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method acting on 
individual unbinned events. 

20 34. A method according to any of claims 31-33 wherein the plurality of detectors consists 
of two such detectors. 

35. A method according to any of claims 31-34 wherein the images are three dimensional 
images. 
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Claims 
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No: 


Claims 




Industrial applicability (lA) 


Yes: 


Claims 


1-35 




No: 


Claims 





2. Citations and explanations 
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claims are fully supported by the description, are made: 

see separate sheet 
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individual, unbinned radiation events which is a concept that does not have anything in 
common with (temporal) binning. The application thus lacks unity. 

V) The concepts of treating individual unbinned events and of treating temporally 
binned events both appear to be novel and inventive. 

VIII) Claims 20 and 21 contain disclaimer features which cast doubt as to the extent of 
protection sought, and thus lack clarity. 

Claims 22 and 23 appear to be lacking essential features, since it is unsufficiently clear 
from the wording of the claims, how exactly the acquired data is utilized in the image 
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CLAIMS 

1. A method of reconstructing tomography images comprising: 

acquiring data on the geometric coordinates of detection of individual radiation events; 
separately distributing a weight of each of the individual radiation events along a line 
5 of night associated v/ith the event detemiined from the acquired data on the geometric 
coordinates of detection of the individual event; and 

iteratively reconstructing the image based on the distributed weights. 

2. A method according to claim 1 wherein the weights are distributed in voxels along the 
10 line of flight and wherein the weight of a particular event is distributed based on the 

probabiHty that an event occurred in particular voxels. 

3. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined based on the position at which the event was detected on a detector and the 

15 acceptance direction of a collimator through which the detector receives radiation associated 
with the events. 

4. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined by the position on a detector on which the event is detected and the location of the 

20 source of radiation associated with the event. 

5. A method according to claim 1 or claim 2 wherein the line of flight associated with an 
event is determined by detection of two coincident photons. 

25 6. A method according to any of the preceding claims wherein iteratively reconstructing 
the image comprises applying an iterative expectation maximization (EM) method on the data 
in sub-sets. 

7. A method according to claim 6 wherein the individual events form the separate sub- 
30 sets. 

8. A method according to claim 6 or claim 7 wherein the sub-sets are formed based on the 
time of acquisition of events. 

35 9. A method according to claim 6 wherein the sub-sets are formed from unrelated events. 
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10. A method of reconstiiicting tomography images comprising: 

acquiring data on the geometric coordinates of detection of individual radiation events; 

and 

^ applying an iterative expectation maximization (EM) method on the data in sub-sets 

which are foimed based on the time of acquisition of the data on the geometric coordinates of 
detection of the events. 

11. A method according to any of claims 6-10 wherein the subsets consist of data having 
10 less than a 180 degree view angle. 

12. A method according to any of claims 6-11 wherein iterations of the EM method are 
performed prior to the acquisition of data having a 1 80 degree angle of view. 

15 13. A method according to any of claims 6-12 wherein iterations are commenced on 
receipt of the first detected event. 

14. A method according to any of claims 6-13 comprising displaying an evolving image 
based on successive iterations iterative method on a display device. 

20 

15. A method according to any of claims 6-14 and including determining if a study should 
be terminated based on the image quality of an image after an iteration. 

16. A method according to any of claims 6-15 wherein intermediate images are filtered 
25 with a smoothing filter between iterations of the EM method. 

17. A method according to any of claims 6-15 wherein intermediate images are filtered 
with a noise reducing filter between iterations of the EM method. 

30 IS. A method according to any of claims 6-17 wherein data is reused in subsequent 
iterations of the EM algorithm. 

19. A method according to any of the preceding claims wherein the image is a three 
dimensional image. 

35 
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20. A method according to any of the preceding claims wherein the iterative method 
comprises reconstmcting from the events without forming two dimensional data sets. 

21. A method according to any of the preceding claims wherein the iterative method 
comprises reconstmcting from the events without forming sinograms for slices of the three 
dimensional image. 

22. A method of reconstructing tomography images comprising: 

acquiring data on the geometric coordinates of detection of individual radiation events; 

and 

iteratively reconstructing a tliree-dimensional image from the unbinned individual 
radiation events. 

23. A method according to claim 22 wherein reconstructing the image comprises utilizing 
an expectation maximization (EM) method acting on individual unbinned events. 

24. A method according to any of the preceding claims wherein the radiation events are 
nuclear emission events and the images are emission tomography images. 

25. A method according to any claims 1-24 wherein the radiation events are positron decay 
events and wherein the images are PET images. 

26. A method according to any of claims 1-24 wherein the radiation events are represented 
by photons which have passed through a subject and wherein the images are transmission 
tomography images. 

27. A method according to claim 26 wherein the radiation events are nuclear 
disintegrations and wherein the images are nuclear transmission tomographic images. 

28. A method according to claim 26 wherein the radiation events are X-rays and wherein 
the images are X-ray CT images. 

29 . A method according to any of the preceding claims wherein the line of flight associated 
with the radiation events form a fan beam. 
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30. A method according to any of claims 1-28 wherein the lines of flight associated with 
the events form a cone beam. 

31. A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on the geometric coordinates of detection of individual positron 

emission tomography events utilizing a plurality of spatially continuous area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method acting on 
individual unbinned events. 

32. A method according to claim 30 wherein the spatially continuous detectors are 
substantially planar detectors. 

33. A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on the geometric coordinates of detection of individual positron 

emission tomography events utilizing a plurality of substantially planar area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method acting on 
individual unbinned events. 

34. A method according to any of claims 31-33 wherein the plurality of detectors consists 
of two such detectors, 

35. A method according to any of claims 31-34 wherein the images are three dimensional 
images. 
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1 . This report has been drawn on the basis of {substitute sheets which have been furnished to the receiving Office in 
response to an invitation under Article 14 are referred to in this report as "originally filed" and are not annexed to 
the report since they do not contain amendments.): 

Description, pages: 

1 -21 as originally filed 

Claims, No.: 
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Drawings, sheets: 

1 -3 as originally filed 
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see separate sheet 

4. Consequently, the following parts of the international application were the subject of international preliminary 
examination in establishing this report: 
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□ the parts relating to claims Nos. , 

V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial 
applicability; citations and explanations supporting such statement 

1. Statement 



Novelty (N) 


Yes: 


Claims 


1-35 




No: 


Claims 




Inventive step (IS) 


Yes: 


Claims 


1-35 




No: 


Claims 




Industrial applicability (lA) 


Yes: 


Claims 


1-35 




No: 


Claims 





2. Citations and explanations 
see separate sheet 

VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the 
claims are fully supported by the description, are made: 

see separate sheet 
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IV) 

Claim 10 and the claims depending thereupon are concerned with a form temporal 
binning (i.e. not spatial binning). The remaining claims are essentially concerned with 
individual, unbinned radiation events which is a concept that does not have anything in 
common with (temporal) binning. The application thus lacks unity. 



V) The concepts of treating individual unbinned events and of treating temporally 
binned events both appear to be novel and inventive. 

VIII) Claims 20 and 21 contain disclaimer features which cast doubt as to the extent of 
protection sought, and thus lack clarity. 

Claims 22 and 23 appear to be lacking essential features, since it is unsufficiently clear 
from the wording of the claims, how exactly the acquired data is utilized in the image 
reconstruction. 
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CLAIMS 

1 . A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; 

distributing a weight of the individual radiation events along a line of flight associated 
5 with the event determined from the acquired data; and 

iteratively reconstructing the image based on the individually reprojected data, 

2. A method according to claim I wherein the weights are distributed in voxels along the 
line of flight and wherein the weight of a particular event is distributed based on the 

10 probability that an event occurred in particular voxels, 

3. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined based on the position at which the event was detected on a detector and the 
acceptance direction of a collimator through which the detector receives radiation associated 

15 with the events. 

4. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined by the position on a detector on which the event is detected and the location of the 
source of radiation associated with the event. 

20 

5. A method according to claim 1 or claim 2 wherein the line of flight associated with an 
event is determined by detection of two coincident photons. 

6. A method according to any of the preceding claims wherein iteratively reconstructing 
25 the image comprises applying an iterative expectation maximization (EM) method on the data 

in sub-sets. 

7. A method according to claim 6 wherein the individual events form the separate sub- 
sets. 

30 

8. A method according to claim 6 or claim 7 wherein the sub-sets are formed based on the 
time of acquisition of events. 

9. A method according to claim 6 wherein the sub-sets are formed from unrelated events. 

35 
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10. A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

applying an iterative expectation maximization (EM) method on the data in sub-sets 
which are formed based on the time of acquisition of the events. 



11. A method according to any of claims 6-10 wherein the subsets consist of data having 
less than a 180 degree view angle. 

12. A method according to any of claims 6-11 wherein iterations of the EM method are 
10 performed prior to the acquisition of data having a 180 degree angle of view. 

13. A method according to any of claims 6-12 wherein iterations are commenced on 
receipt of the first detected event. 

15 14. A method according to any of claims 6-13 comprising displaying an evolving image 
based on successive iterations iterative method on a display device. 

15. A method according to any of claims 6-14 and including determining if a study should 
be terminated based on the image quality of an image after an iteration. 

20 

16. A method according to any of claims 6-15 wherein intermediate images are filtered 
with a smoothing filter between iterations of the EM method. 

17. A method according to any of claims 6-15 wherein intermediate images are filtered 
25 with a noise reducing filter between iterations of the EM method. 

18. A method according to any of claims 6-17 wherein data is reused in subsequent 
iterations of the EM algorithm. 

30 19. A method according to any of the preceding claims wherein the image is a three 
dimensional image. 

20. A method of reconstructing tomography images comprising: 



5 



acquiring data on individual radiation events; and 
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iteratively reconstructing a three-dimensional image from said individual radiation 
events without producing a stack of two dimensional data sets. 

21. A method of reconstructing tomography images comprising: 
^ acquiring data on individual radiation events; and 

iteratively reconstmcting a three-dimensional image without producing individual 
sinograms for slices of the three dimensional image. 

22. A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; and 

iteratively reconstructing a three-dimensional image utilizing the individual radiation 
events without spatially binning of the events. 

23. A method according to any of claims 20-22 wherein reconstructing the image 
comprises utilizing an expectation maximization (EM) method acting on individual unbinned 
events. 

24. A method according to any of the preceding claims wherein the radiation events are 
nuclear emission events and the images are emission tomography images. 

25. A method according to any claims 1-24 wherein the radiation events are positron decay 
events and wherein the images are PET images. 

26. A method according to any of claims 1-24 wherein the radiation events are represented 
by photons which have passed through a subject and wherein the images are transmission 
tomography images. 

27. A method according to claim 26 wherein the radiation events are nuclear 
disintegrations and wherein the images are nuclear transmission tomographic images. 

28. A method according to claim 26 wherein the radiation events are X-rays and wherein 
the images are X-ray CT images. 

29. A method according to any of the preceding claims wherein the line of flight associated 
with the radiation events form a fan beam. 
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30. A method according to any of claims 1-28 wherein the Unes of flight associated with 
the events form a cone beam. 

31. A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on individual positron emission tomography events utilizing a plurality 

of spatially continuous area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method. 

32. A method according to claim 30 wherein the spatially continuous detectors are 
substantially planar detectors. 

33. A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on individual positron emission tomography events utilizing a plurality 

of substantially planar area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method. 

34. A method according to any of claims 31-33 wherein the plurality of detectors consists 
of two such detectors. 

35. A method according to any of claims 31-34 wherein the images are three dimensional 
images. 
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AMENDED CLAIMS 

[received by the International Bureau on 12 Oune 1998 (12.06.98); 
original claims 1, 31 and 33 amended ; 
remaining claims unchanged (2 pages)] 

1 . A method of reconstructing tomography images comprising: 
acquiring data on individual radiation events; 

separately distributing a weight of each of the individual radiation events along a line 
of flight associated with the event determined from the acquired data; and 
iteratively reconstructing the image based on the distributed events. 



2. A method according to claim 1 wherein the weights are distributed in voxels along the 
line of flight and wherein the weight of a particular event is distributed based on the 

10 probability that an event occurred in particular voxels. 

3. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined based on the position at which the event was detected on a detector and the 
acceptance direction of a collimator through which the detector receives radiation associated 

15 with the events. 



4. A method according to claim 1 or claim 2 wherein the line of flight of an event is 
determined by the position on a detector on which the event is detected and the location of the 
source of radiation associated with the event. 

20 

5. A method according to claim 1 or claim 2 wherein the line of flight associated with an 
event is determined by detection of two coincident photons. 

6. A method according to any of the preceding claims wherein iteratively reconstructing 
25 the image comprises applying an iterative expectation maximization (EM) method on the data 

in sub-sets. 



7. A method according to claim 6 wherein the individual events form the separate sub- 
sets. 

30 

8. A method according to claim 6 or claim 7 wherein the sub-sets are fomied based on the 
time of acquisition of events. 

9. A method according to claim 6 wherein the sub-sets are formed from unrelated events. 

35 
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30. A method according to any of claims 1-28 wherein the lines of flight associated with 
the events form a cone beam. 

5 31, A method of reconstructing positron emission tomography (PET) images comprising: 
acquiring data on individual positron emission tomography events utilizing a plurality 

of spatially continuous area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method acting on 

individual unbinned events. 

10 

32. A method according to claim 30 wherein the spatially continuous detectors are 
substantially planar detectors. 

33. A method of reconstructing positron emission tomography (PET) images comprising: 
15 acquiring data on individual positron emission tomography events utilizing a plurality 

of substantially planar area detectors; and 

reconstructing the image utilizing an expectation maximization (EM) method acting on 
individual unbinned events. 

20 34. A method according to any of claims 31-33 wherein the plurality of detectors consists 
of two such detectors. 

35. A method according to any of claims 31-34 wherein the images are three dimensional 
images. 
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(i) 


the claim 


is unchanged; 


(K) 


the claim 


is cancelled; 


(Hi) 


the claim 


is new; 


(iv) 


the claim 


replaces one or more daims as filed; 


(v) 


the claim 


is the result of the division of a claim as filed. 
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